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Orbital elementary excitations as probes of entanglement and quantum phase 
transitions of collective spins in an entangled Bose-Einstein condensate 
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A mixture of two species of pseudospin-i Bose gases exhibits interesting interplay between spin 
and orbital degrees of freedom. Expectation values of various quantities of the collective spins of the 
two species play crucial roles in the Gross-Pitaevskii-like equations governing the four orbital wave 
functions in which Bose-Einstein condensation occurs. Consequently, the elementary excitations of 
these orbital wave functions reflect properties of the collective spins. When the coupling between 
the two collective spins is isotropic, the energy gap of the gapped orbital excitation peaks, while 
there is a quantum phase transition in the ground state of the effective Hamiltonian of the two 
collective spins, which have previously been found to be maximally entangled. 
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I. INTRODUCTION 

A current trend in condensed matter physics is the 
interplay between spin and orbital degrees of freedom. 
It is interesting to explore this topic in the realms of 
Bose-Einstein condensation (BEC), as physical effects on 
single-particle levels are often enhanced by Bose statis- 
tics. A possible avenue is spinor Bose gases [H-Q- The 
perspective is more broadened if we consider a mixture of 
two distinct species of spinor Bose gases, where there are 
several coupling strengths for the collective spins of the 
two species, which depend upon the orbital wave func- 
tions. This situation leads to coupling between spin and 
orbital degrees of freedom. For simplicity, let us focus 
on a mixture of two species of pscudospin-i Bose gases 
with both intraspecies and interspecies spin-exchange in- 
teractions [B-l^ ■ Interspecies spin-exchange causes inter- 
species entanglement and the so-called entangled Bose- 
Einstein condensation (EBEC) Q, which is different from 
a mixture of two species of spinless atoms, where the 
two species are disentangled and individually undergoes 
BEC. EBEC amplifies quantum entanglement from indi- 
vidual particles to macroscopic condensates. This ground 
state bears some analogies with the lowest energy state of 
a single species of pseudospin-i atoms in a double well @ 
or occupying two orbital modes d, [13, [HI , but there are 
also differences due to the fact that two atoms of dis- 
tinct species are distinguishable and thus the numbers of 
atoms are respectively conserved. 

It has been shown that the interspecies entanglement 
in the ground state is maximal at the isotropic parame- 
ter point of the effective Heisenberg coupling of the two 
collective spins, and that the larger the particle numbers 
of the two species, the steeper the entanglement peak [^. 
In this brief report, after elucidating that the entangle- 



ment peak is indeed located at a quantum phase transi- 
tion point, we study how the ground state of the collec- 
tive spins affects the elementary excitations of the orbital 
wave functions in which EBEC occurs. Furthermore, we 
find that in the vicinity of this quantum phase transition, 
the energy gap of the gapped orbital elementary exci- 
tation is strikingly different from those of disentangled 
ground states. Away from the quantum phase transition 
point, the elementary excitations tend to approach those 
of a disentangled ground state. 



II. THE MODEL 

Consider a dilute gas of two distinct species of atoms. 
Each atom has an internal degree of freedom represented 
as a pseudospin-i, with z-component basis states f and 
|. For a single species of pseudospin-i gas, it was ar- 
gued that the conservation of a small total spin in the 
cooling process invalidates the single orbital mode ap- 
proximation [lol . [Tl| . In our system, in contrast, un- 
der the single orbital mode approximation, the total spin 
of the mixture can still be arbitrarily small due to the 
distinguishability of the two species. Therefore, there 
is no reason against the usual single orbital mode ap- 
proximation, which works well in most of the BEC sys- 
tems. Moreover, a field theoretical approach without us- 
ing single orbital mode approximation confirms its valida- 
tion 



12|. Therefore, for each atom of species a{a = a, b) 
and pseudospin a {a —'[,],), we can safely assume that 
only the single-particle orbital ground state (jiaa (r) is oc- 
cupied. Then the many-body Hamiltonian can be writ- 



ten as H = E a, a U.N,, + 5 

Kf^,Na,Nb,'+Keia\a^blb^ + ala^b\b^), whe 
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^aa' aaJ-^ba' + ^^e[a\ait}'^t}-^ + a'^afb:^bi), Where aa 
denoted annihilation operator associated with 4)a,(r) 
of species a, N^a = cel-a, is the corresponding 
number of atoms For each specie a, the to- 



tal particle number 



N, 



at 



Nai is a constant. 
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The coefficients K' s are shorthand for K'^^^^^^^^ 



9\ 



(Ti (72 CTs (74 



i4>aai<P'r, where 



ja\a2 <73C74 



27i-/i2C?f^,^3^,//i„/3, with CfcT2<T3<T4 being the scattering 
length for the scattering in which an a-atom flips its pseu- 
dospin from to cri while an /3-atoni flips its pseudospin 
from (T3 to (72, pLap = 'm,a'mp/{ma + mp) is the reduced 
mass. For intraspecies scattering, K^^ = K^"^^,K^^ = 
3 for a ^ a. For interspecies scattering, 
= K°-^ — K°-^ Pnli f ^ f — 
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K^^I2, where 
is the single-particle energy. 

For simplicity, we assume that the scattering lengths 
satisfy g'^'^^^ = g^^^^ = g^, gf„„„ = gs, gfgg^ = ga, 



g'^l^^ = = 9e- Moreover, we focus on the uniform 

case (j)aa = where is the volume of the system. 

The total spin operator of species a is Sq, = a\scrcria„i , 
where s^a' is the single spin operator. In terms of S^, the 
Hamiltonian can be transformed into that of two coupled 
giant spins Sa = Na/2 and Sb = Nb/2, 
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FIG. 1: Fidelity susceptibility x as a function of 77 e 
Qd) for A'^a = Nb = N . The left inset shows how r]. 
with A'' and the right inset shows how xiv^ 



■max varies 
varies with A''. 



H — J±{SaxSbx + SaySby) + JzSazSbz + Eq (1) 



where = = ^ 

2S1 



-ad) 



Eo 



Na{Ng 



-9a 



2S1 



-9b 



'-{9s 



gd) is a constant. Without loss 



of generality, let Sa> Sb- We focus on antifcrromagnetic 
couplings 5e > and gs — gd > 0- 



III. QUANTUM PHASE TRANSITION 

There is a quantum phase transition at the parameter 
point ge = 9s — 9d- The ground states are qualitatively 
different in the limits of ge 3> gs — 9d and g,, <^ ga — gd- 

In a mean field approximation, the ground state is dis- 
entangled between the two species, and can be written 
as \G)mf = 
(e-^^o/^cos^ 
field energy E^p 

^-^^^{Saz){Sbz). 



n 



I t)fc + e"^''/^sin|i| \.)b)^\ with mean- 

. ^f{{Sax){Sb.) + {Say){Sby))) + 

where a constant is neglected. 



(Sax) = ^sin6'aCOSV3Q, (Say) = ^sine'^sinvsa, {Sa 

N, 



-^cos9a, (a 



a,b). 

For ge < gs — gd, Emf is minimal when 0a = while 
9b ~ TT or = TT while 9b = 0, that is, the two spins are 
antiparallel and are along z axis. 

For ge > gs — gd, Emf is minimal when 9a = 9b = tt/2 
while ipb — (pa + TT, that is, the two spins are antiparallel 
and are on x — y plane. 

One can also use the so-called fidelity susceptibility to 
analyze the QPT [H, For a Hamihonian n{7]) = 

"Ho + rfJii, where 77 is a driving parameter, the ground 
state fidelity is defined as F = |(i/'o('? + <5?7)|'/'o('7))|- For a 
nondegenerate ground state, the fidelity susceptibility is 

XVI) — mn5^_>o -3;^ - l^n^o [E„in)-Eo{vW ' ^"'^^'^ 
|'0o(77)) and |V'ra('7)) are the ground and excited states of 
T-L^q), respectively. In our Hamiltonian, "Hq = JzSazSbz, 

T-Ll = Jz{SaxSbx + SaySby), rj = Jj_/Jz = ge/ [gs ~ 9d)- 



We calculate the fidelity susceptibility for the case of 
Na = Nb, in which the ground state is non-degenerate. 
In FigUJ i]max is the position where the fidelity suscepti- 
bility is maximal, r^rnax ^ l~5.33i9N-^-^^^^, the fidelity 
susceptibility becomes sharper in the vicinity of r]max 
when the number of atoms increases, x ~ TyS-S^i j^^i^ 
appearance of the peak implies profound change of the 
ground state at rjmax- In the thermodynamic limit, rjmax 
approaches the maximum 1, where QPT takes place. 

The result of the fidelity susceptibility is very much 
consistent with the previous result of the entanglement 
entropy Q. First, the maxima of both quantities are lo- 
cated at 7/ = 1. Second, the trends in approaching of 
the maxima are also consistent. The larger the parti- 
cle number N , the quicker both quantities approach the 
maxima. Third, both quantities display asymmetries be- 
tween 77 > 1 and 77 < 1, with the approaches of both 
quantities to maxima being quicker on the side of 77 < 1 . 



IV. ELEMENTARY EXCITATIONS OF THE 
ORBITAL WAVE FUNCTIONS 



Now we show how the spectra of the elementary exci- 
tations of the orbital wave functions vary with the many- 
body ground state. These single-particle wave func- 
tions are governed by the the generalized time-dependent 
Gross-Pitaevskii-like equations, which are determined by 
minimizing the energy, i.e. the expectation of the Hamil- 
tonian under the many-body ground state, as a func- 
tional of the orbital wave functions [H . Previous calcula- 
tions have been focused on the many-body singlet ground 
state 0], here we make a more general consideration. 

From the Hamiltonian, wc obtain the Gross-Pitaevskii 
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equations 

+ {Nl^ - iV„^)ga|0aaP + (iVa<TVaff)5a|0a<Tp (2) 
+ {NcaNp„)gs\(t)[3a? + {Nc.aNps)gd\(t)ps?}Ka 

where /? 7^ a, ct 7^ cr, (O) represents the expectation value 
of the operator O in the many-body ground state. 

The particle number operators whose expectation val- 
ues appear in the Gross-Pitacvskii-likc equations can be 
represented in terms of collective spins of the two species, 
because of the relations Nacr = + i]aSaz, with 

77^ = 1 and rii = -1, and {a'^asPlPa) = {SaxSbx + 
SaySby). Therefore, (iV„,) = + 7j,{Saz) , {Nl„) = 

Nil A + 77,iV„(5„,) + {Si,), {NaaN^a) = Nl/A - {Si,), 

{NaaNba) = {NaNb)/A + {SazSbz) + + 

?7,(iV6/2)(5a.), = {NaNb)lA - {S^,Sp,) + 

ils{Na/2){S,3,) + VaiNp/2){Saz)- In the special case of 
many-body singlet ground state, we have Na = Nb = N 
and {Saz) = {Sbz) = {Sz) = 0, leading the special form 
of Gross-Pitaevskii-like equation considered previously. 

With these spin quantities as the coefficients, the 
Gross-Pitaevskii-like equations reflect coupling between 
spin and orbital degrees of freedom. As the ground state 
varies with the parameters, so are the expectation values 
of these spin quantities. 

For simplicity, we focus on the uniform case in ab- 
sence of external field U = 0, the lowest energy wave 

n — "i- l-i- a cT t / ^ 

functions are thus (f>^^ = - — -= — , where is the 



where p = N/n, A, = -{SazSbz)/{Ny4), = 
-{SaxSbx + SaySby) /{N'^l 4), and we have set {Nl^ - 
Nc.a)l{N^a) « {Nl,)l{N^,). 

It can be found that there are four energy spectra of 
elementary excitations, 

fiwi,2 = -^^V,TVWif+^2. (7) 

^^M-7f\/r3Tv/(r3)2 + r4, (8) 



chemical potential determined by substituting (jr^^ in 
@, (A^a<.)/Xa. = {Nl, - A^a.)t + {N^aNo.^)a^ + 
{Nc..Np„)% + {N^,Np^)^ + (4a<j4/3,)ff . 
Now we consider elementary excitation, i.e. deviation 

(j^cta- Then the Gross- Petaevskii-like equation yields the 
Bogoliubov-de Gennes-like equations 

{N^„){hu - E'^)u^„ = 

{NI„ - Naa)^{Uaa ~ Vaa) + {NaaNaa)^{Uaa ~ Vaa) 
+ {NaaNbcr)^{upa ~ VfSa) + {NaaNpa)^{uf!a ~ Vpg) 
+ {alagl3ll3a)^{ui:Sa + UaS - Vpg ~ Uaa), 

(3) 

{N^^){-hw-E^)v,,,, = 

{NI^ - Naa)^{Vaa ^ "acr) + {NaaNaa)^{Vaa ~ Uaa) 
+ {NacrNbcr)^iv/3cr ~ Up^r) + {NaaN i3a)^{Vf3iT " Ujia) 
+ {a\aafilfia)^{vi3a +Vas ~ Upg - Vaa)- 

(4) 

^herei?," = £- 

For simplicity, in the following we focus on the ground 
states with Na = Nb = N, total 5*2 = 0, and thus {Saz) = 
0, {SD = -{SazSbz). Therefore {N^a) = N/2, {N^^) = 

{NaaNpg) = N^A- {SazSbz), {NaaNag) = {NaaNba) = 

N'^/4:+ {SazSbz)- Hence © and (g]) can be written as 



(5) 



(6) 

I 

{E^+2pga){E'g + 2pgb)}, = E^[E^ + 2p{gaXz+geX±)] + 
E^ [E^ + 2p{gbXz +5e A_L )] + 2p^g^X^ [{ga+9b + gs+gd)Xz + 
ge\^+gd~gsl = A[E^^El + pg,\A.{E''q+ El)]{p'[gs- 

9d - {9s + gd)Xz?- [E^ + p{2gaXz + 5e A± )] [E^^ +p{2gbXz + 
<?eAj_)]}. If ge = 0, the Hamiltonian would have four U{\) 
symmetries which result in four gapless energy spectra. 
Now that ge 7^ 0, there are only three U{\) symmetries, 
three energy spectra are gapless, and the fourth energy 
spectra opens a gap 0]. Indeed, when g = and ge 7^ 0, 

^1,2,3 = 0, hijji ^ 0. 



J 



{hhJ - E^)Uaa \pga{l + X,){Uaa - Vaa) + ^P9a{i " X,){Uag - Vag) 

+ ^pgs{l - \z){upa - Vp„) + \p9d{'^ + Xz){upg ~ Vpg) - \pge\±{upa + UaS - Vpg - Uaa), 
{-huj - E^)Vaa = \p9a{^ + K){Vaa " Uaa) + " X,){Vag - Uag) 

+ \pgs{l - \z){vpa - Upa) + ^P9di^ + K){vpg - Upg) - ^pge\±{vpa + Vaa " Upg - Vaa) = 0, 



where — 's are for hjji and hw^, respectively, -l-'s are for 
hu}2 and hhj^, respectively, Fi = E^{E^ + 2pga) + E^^{E^q + 
2pgb), F2 = AE^E^^ip^gs + gd ~ (ff. - gd)Xz - ffeA^l^ - 



We now study the behavior of the energy gap A of the 
forth elementary excitation spectrum as a function of ge 
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and gs — g^. It is easy to write down 

A= ^ yrs + y/TfTTI 

= P V2geA_L [{9a + 9b + gs + 9d)^z + ,9eA_L + gd - 9s] 

(9) 

Obviously A 7^ only '\i 9e 7^ 0. Hence interspecies 
spin exchange is necessary for opening a gap in an orbital 
excitation. 

We have numerical calculated \z and Aj^ for various 
values of ge and gs — gd, and then obtained A by using 
(|9]). It is found that when N 00, 

ffe < 9s — gd ■ -^2 ^ 1, A_L 0, hence A — >• 0, 

9e = 9s - 9d ■■ K 1/3, X± 2/3, (10) 

9e> 9s- 9d- ^ 0, A_L 1- 

It can be seen that the numerical result at g^ = 9s — gd is 
the same as the exact result. A_l — > implies that the ef- 
fect of spin exchange diminishes, consequently wa — >■ wi, 
a;4 — > aj2, the excitations become two doubly degenerate 
ones. Note that this is under the condition that each 
scattering strength is symmetric for the two pseudospin 
states of each species. 

As shown in Fig. ^ the larger the numbers of atoms 
of the two species, the closer is the gap to that under the 
disentangled ansatz, except in the vicinity the critical 
point, where the gap varies rapidly with ge/igs — gd)- 
At critical point, as the atom numbers increase, the gap 
quickly saturates to a non-zero value. 

Under a disentangled mean-field ground state, one has 
= 1 and Aj^ = for g^ < g^ — gd, and Xz ~ and 
A_L = 1 for ge > gs — gd- Therefore, far away from the 
critical point ge — gs ~ gd, Xz and A_l, and thus the 
actual elementary excitations, are close to those under 
the disentangled mean-field states. But the disentangled 
ansatz clearly fails in the vicinity of the critical point 
5e — gs — gd, as it tells that Xz and Aj^ are arbitrary 
non- negative values satisfying A + Aj^ = 1. 

V. SUMMARY 

A mixture of two pseudospin- i Bose gases with inter- 
species spin exchange displays interesting interplay be- 
tween spin and orbital degrees of freedom. The many- 
body Hamiltonian is simplified to an anisotropic Heisen- 
berg coupling between the two collective spins of the 
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FIG. 2: The energy gap, scaled by \/2p{ga — Qd), as a function 
of ge/{ga — gd), for different values of N. Points represent 
numerical results, while the solid line is from the mean-field 
disentangled ground state. The parameter values used satisfy 

{ga + gb + gs+ gd)/{gs - gd) = 5.0 



is critical point of quantum phase transition, where the 
interspecies entanglement and fidelity susceptibility also 
peak. These properties should be carried over to a spin-1 
mixture [l6j . 
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